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Abstract 

Recent developments in the Symanzik improvement program for lattice 
QCD are reviewed. 

O , 

^ ■ 1 Introduction 

. In this lecture I review recent trends in the effort to develop highly improved 

discretizations of QCD for use on coarse lattices. I focus here on the perturbative 
Symanzik approach to this problem; other approaches are discussed by other 
lecturers at this meeting. The ultimate goal of all this activity is too allow us 
to simulate with larger lattice spacings, at greatly reduced cost. The savings 
that result from larger lattice spacings can be used to extend the reach of lattice 
QCD to far more difficult problems, well suited to our largest supercomputers. 
At the same time, there is a wealth of simpler simulation problems that are now 
. accessible to anyone with a modern workstation or personal computer. 

$H ' It is somewhat artificial to isolate a single approach like Symanzik improve- 

ment. Symanzik improvement is but one of a large set of tools — Symanzik and 
on-shell improvement, perturbation theory and tadpole improvement, blocked 
fields and perfect actions, nonperturbative tuning and MCRG, and so on — that 
we can use to design improved discretizations. These tools are closely related, 
and not mutually exclusive. Solving QCD is much too difficult for us to ignore 
any particular set of tools; we need every trick available. 

I begin, in Section ^ with a review of perturbative improvement and tadpole 
improvement. I discuss variations on the canonical approach, including non- 
perturbative tuning, and I discuss what we now know about the size of finite-a 
errors. Weak-coupling perturbation theory plays a central role in all practical 
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approaches to improved discretization. I discuss the status of lattice perturba- 
tion theory, as well as how to implement and test it, in Section [| In Section |] 
I discuss specific, recent developments concerning various actions. Finally I 
summarize in Section H. 



2 Perturbative and Tadpole Improvement 

2.1 The Method 

The procedure for perturbatively improving a lattice QCD operator or action 
involves three steps jjj : 

1. Classical Improvement: The classical (tree- level) lattice action or opera- 
tor is corrected through the desired order in the lattice spacing a. For 
example, the lattice gluon action 

s = -p {ci p »A x ) 

X,fl>lS 

+c 2 (RwA x ) + R ™A X ))} (!) 



equals 

|d 4 *±TrF^) + 0(a 4 ) (2) 

for classical fields when 

/3 = 6/.g 2 , ci = 5/3, c 2 = -l/12. (3) 



Here P^y is one third the real part of the trace of the lxl Wilson loop 
in the \iv plane, while is the same for a 2 x 1 loop with two steps 

in the ^ direction and one in the v direction. We derive this result by 
substituting the definition of the lattice link, in terms of the continuum 
vector potential, 



Un(x) = Pexp 



(4) 



into the lattice action, and expanding in powers of the lattice spacing a. 

2. Tadpole Improvement: In the quantum theory, the couplings are renor- 
malized by quantum fluctuations, and new operators are required. The 
dominant renormalizations are due to tadpole diagrams caused by the 
(nonlinear) gluon link operator U^. These lattice artifacts are largely can- 
celled by tadpole improvement: every link in the classical lattice theory is 
divided by the mean link |l[] : 

U,{x) (5) 
u 
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The mean link uo is a number between and 1; it is computed self- 
consistently, using a Monte Carlo simulation, from the mean value of 
■^TiU^. This last quantity is gauge dependent and can be made arbitrarily 
small by changing the gauge. Thus we choose the gauge that maximizes 
(Tr[/ M ), which is Landau gauge: 

u = (0|iTr[/ M |0) L G- (6) 

This gauge minimizes the tadpoles in uq, pushing it as close to 1 as pos- 
sible; any tadpole contribution that remains cannot be a gauge artifact. 
Our gluon action, Eq. (|l|), is tadpole improved by setting 

The tadpole correction here is substantial at lattice spacings of .25-.4fm 
since uo is .7-. 8 and is raised to large powers. 

At small lattice spacings uo is well approximated by the fourth root u p of 
the plaquette operator, which is easier to compute than the mean link in 
Landau gauge: 

u ~u p = (OIP^IO) 1 / 4 . (8) 

3. Radiative Corrections: A tadpole-improved classical action or operator 
is an efficient starting point for a systematic weak-coupling expansion. 
Each coupling constant has an expansion in powers of the QCD coupling 
constant ay- 

(H = -jjt (cj +c,i a v (qf) H ) (9) 

V 

The scale q* is typically between 0.5/a and 2/a, and is easily computed 
using one- loop perturbation theory [^| . The value of the coupling con- 
stant is also easily measured in the simulation |2j. New operators, beyond 
those necessary in the classical analysis, are usually needed when radia- 
tive corrections are included. The perturbative expansions of the coupling 
constants are computed by matching: low-energy physical quantities, scat- 
tering amplitudes for example, are computed perturbatively in both the 
lattice theory and the continuum, and lattice couplings are tuned so that 
the two calculations agree order-by-order in the coupling (to a given order 
in the lattice spacing). The perturbation expansion is usually much more 
convergent with the tadpole factor, l/uVj*, removed as above. Indeed, for 
many applications, only the tree-level term in the expansion is necessary 
after tadpole improvement. 

The procedure just outlined is a double expansion in powers of the lattice 
spacing and of the QCD coupling constant. Carried to all orders it becomes 
infinitely complicated; in practice, however, only a few low orders of correction 
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are needed to reduce errors to a few percent. Recently we have started to obtain 
high-quality data on the size of finite-a errors. An example is the SCRI analysis 
of errors in the hadron spectrum for Wilson and SW quarks || . By going to 
large a, they made the errors large and unambiguous; they tested the algorithms 
"to failure." They find that the relative errors are roughly (a/ao) n with n=l 
for Wilson quarks and n = 2 for SW quarks, and a w 1-2 fin. The powers are 
as expected, and so is the scale a for the errors §. 

Such information tells us how much improvement we need at a given lattice 
spacing. At a = .4fm, for example, a/ao ~ -4 and ay ~ -3. Thus, if we 
want errors smaller than 5-10%, we should include in our quark action tree- 
level corrections through order a 3 and radiative corrections through order a ay 
(that is, radiative corrections for only the order a term, and then only through 
one-loop). At a = .lfm, on the other hand, only the order a and, perhaps, the 
order a ay corrections are necessary for 1% accuracy (in addition to 4 6 times 
as much computing). Note that as a increases higher powers of a become more 
important faster than higher powers of ay . These arguments are approximate, 
but they provide a starting point for a realistic assessment of what is needed. 

Radiative corrections (beyond tadpole improvement) are unimportant for 
most terms in quark and gluon actions. This is because the radiative corrections 
are in terms already suppressed by some power of the lattice spacing; radiative 
corrections for terms of order a are generally absorbed into masses or other 
parameters that are tuned numerically. The situation is different, however, 
for lattice operators like currents. Typically the desired continuum current is 
obtained from a lattice current using 

Jco„t = Z{a v ) J ut + ll (a v )a5J 1 + 0{a 2 ). (10) 

While radiative corrections may be unimportant in the 0(a, a 2 . . .) corrections, 
the Z-factor appears in order a° and therefore must usually be determined 
through one or two orders in ay to achieve useful accuracy. The problem is 
even more accute for operators that mix through radiative corrections with 
operators of lower dimension, because of 1/a divergences; this can be alleviated 
to some extent by working with larger a's, where 1/a is then smaller. 

An alternative to tadpole improvement worth investigating is to replace the 
link operators in quark actions and currents by smeared links: for example, 

U» (11) 

where A^ 2 ) is the discretized gauge-covariant laplacian, n is a small integer, and 
e a smearing parameter. This introduces a factor like exp(— e a 2 q 2 ) into the gluon 
propagator which reduces the effective ultraviolet cutoff for gluons coupled to the 
quarks. Since tadpoles are quadratically divergent, even a modest reduction in 
this cutoff greatly reduces their contribution and makes tadpole improvement 
unnecessary. Substituting Un for in a lattice operator introduces a new 
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a 2 error; this can be removed by including an additional factor of (1 — ea 2 A^') 
in the definition of U^. 

2.2 Nonperturbative Tuning 

A variation on the standard procedure is to tune some of the couplings nonpcr- 
turbatively, using simulations. There are two standard approaches. One is to 
tune a coupling constant to restore a continuum symmetry. In many cases a 
leading correction term breaks a continuum symmetry, like Lorentz invariancc, 
and is the only term of its order (in a) to do so. Then the coefficient of that 
term can be determined by adjusting it until the corresponding symmetry is 
restored in the simulation results. For example, there is only one independent 
dimension-six operator that breaks rotation invariance and contributes to the 
0(a 2 ) corrections in a gluon action. Thus we can adjust coupling C2 in our 
lattice gluon action, Eq. (^), until, for example, the simulation gives a static 
quark potential with V(2r,2r,r) = V(3r, 0,0) for some r. The tuned coupling 
then automatically includes perturbative corrections to all orders in ay, as well 
as any nonperturbative contributions!!! The tuned coupling is not exact, how- 
ever, since the quantities used to tune it have finite-a errors caused by the 
action. In our example, the gluon action is accurate only through order a 2 . 
Thus V(2r, 2r, r) — V(3r, 0, 0) has order a 4 errors; tuning it to zero induces or- 
der a 4 errors in C2- Obviously, such truncation errors are smaller for more highly 
corrected lattice theories. Also some care is necessary in choosing the quantitiy 
to be tuned; one wants a quantity for which finite-a errors are minimal or small. 

Other examples of couplings that can be tuned using symmetry restoration 
include: the clover coefficient in the Sheikholeslami-Wohlert (SW) quark action, 
using chiral symmetry [5j; the coefficient of the tpA^ - -ftp correction in the D234 
actions, using Lorentz invariance (see below); and n t /n s in the Fermilab quark 
action for massive quarks, using Lorentz invariance. 

The second approach to nonperturbative tuning is to adjust couplings so 
that physics remains the same at different lattice spacings. Thus, for example, 
one might adjust the clover coefficient in the SW action so that the rho mass 
scales correctly from lattice spacing to lattice spacing. This is harder than 
tuning for a symmetry since simulations at several lattice spacings are required, 
but it is not impossibly difficult if only one or two couplings are needed. One 
loses one "prediction" from the simulation for each coupling that is tuned this 
way. Again the couplings that result include all radiative corrections and any 
nonperturbative effects; and, again, the couplings have finite-a errors caused by 
the action. 

1 Any short-distance quantity could have nonperturbative contributions. Usually these are 
small compared with the perturbative contributions, which is why perturbative QCD works 
so well. Note that lattice coupling constants have renormalon contributions because of finitc-a 
errors in the lattice action; but these contributions are automatically suppressed to the same 
level as other finite-a errors caused by the action, and so can be ignored. 
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Nonperturbative tuning is particularly useful for leading-order corrections, 
like the order-a term in quark actions, where there are only one or two couplings 
to tune. Where it can be used, it is also useful for currents since otherwise 
one-loop and sometimes two-loop radiative corrections must be computed. For 
finite-a corrections beyond leading order, however, tadpole improvement and 
perturbation theory are likely to be the most feasible approach. Fortunately, 
tree-level perturbation theory is probably adequate for most correction terms 
in order a 2 and higher. 



Weak-coupling perturbation theory plays an essential role in all current ap- 
proaches to improved lattice QCD actions and operators. It tells us that im- 
proving the classical theory is a good starting point for the quantum improve- 
ment program. It tells us which corrections are most important; as a result 
we need only a small number of correction terms. In the perturbative improve- 
ment schemes I discuss here, perturbation theory is the tool for calculating the 
coupling constants. But can we trust perturbation theory? 

The short answer is yes. Wherever careful tests have been made, using per- 
turbative expansions through second or third order, for quantities that are easy 
to measure accurately in simulations, perturbation theory has been successful. 
For example, in f| four values for the QCD coupling constant ap are extracted 
using third-order expansions for four different Wilson loops — W\i, W12, W13, 
and W22 — at lattice spacing awO.lfm (/3 = 6). These agree with each other 
to within the estimated truncation errors, 0(a P ) ~ 1%; and they agree to the 
same precision with the value obtained from the (very different) Schrodinger 
Functional approach |8| . This is a high-precision test of QCD perturbation the- 
ory, much more stringent than any test based on continuum quantities. Another 
three-loop test is coupling constant evolution: As illustrated in Fig |], the run- 
ning of the QCD coupling constant is well described by the perturbative beta 
function at least down to scales of order 5 GeV; less accurate tests show this is 
true even down to scales of order 1 GeV Q . 

Usually when we use perturbation theory, we truncate the expansion at 
some finite order in ay. We should always estimate the truncation error that 
results. Sometimes an apparent disagreement between perturbation theory and 
simulation results is negligible when compared with the truncation error, and 
is therefore of little significance. For example, the critical k for SW quarks in 
tadpole- improved perturbation theory is known only through first-order [^[ (t| : 



3 Perturbation Theory 



3.1 Does It Work? 
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lattice spacings, plotted versus the typical momentum, q\s — 3.4/a, in the 
plaquette; lattice spacings are from the 1S-1P upsilon splittin in NRQCD. The 
solid line is predicted by the three-loop beta function. 

where u p is the plaquette mean link. Thus perturbation theory should only 
be correct up to relative errors of order ±\ay. In Figure |5[ we compare the 
first-order result, including error bars for higher-order terms, with nonperturba- 
tive results obtained from simulations Q . The nonperturbative results are more 
accurate, but perturbation theory agrees well — that is, it agrees to within trun- 
cation errors. 

It is equally important, when using nonperturbative tuning, to estimate the 
finite-a errors due to truncation of the lattice action. The clover coefficient c sw 
in the SW action can be tuned nonperturbatively using chiral symmetry J5| . 
The tuned c sw , however, has truncation errors of relative order a/ag for some 
ao around 1-2 fm, at least for larger lattice spacings (at smaller lattice spacings 
one can compute the order a term, using perturbation theory, and remove it). 
The tuned coupling can be compared with the perturbative result || , 



which also has truncation errors, of relative order ±la v . Scale g* w is currently 
unknown, but is likely around .5 /a— 2/a. In Fig. || we compare the nonper- 
turbative and perturbative results, indicating the uncertainties in each, for a 
range of /3's. (For this we guess a = 2fm and q* w = 1 /a). The perturbative 
and nonperturbative c sw 's agree well, and are roughly equal in accuracy; each 
reinforces our confidence in the other. Note that, since c sw is the coefficient of 
an O(a) correction, the different c sw 's at /3 = 6 give identical hadronic physics 
to within 1-2%. 

We also include results in the last figure from first-order perturbation theory, 
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- [1 + 0.20 a v (q* sw ) + 0(a 2 v )] , 



(13) 
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6 6.2 6.4 6.6 6.8 

Figure 2: Critical n for SW quarks from "Monte Carlo simulations and from 
first-order perturbation theory. Statistical errors are negligible for the Monte 
Carlo results; relative errors shown for the perturbation theory are ±la v . 



but with the Landau-gauge mean link uq instead of u p : 

c sw = -L [1 + 0.46 a v + 0(a 2 v )] . (14) 

Uq 

While uq and u p are almost identical through order ay, they differ in higher 
order. Thus differences between the two options for tadpole improvement indi- 
cate the potential size of order a.y corrections; these differences are 5-10 times 
smaller if a v corrections are included in the two c sw expansions. 

Perturbation theory triumphs in these and most other examples; it is usually 
neither more nor less accurate than expected (from perturbation theory). Nev- 
ertheless there will be situations where perturbation theory fails to converge. 
In QED, for example, the decay rate for orthopositronium is proportional to 
a 6 (1 — 10.30/71") 0. This expansion would be useless if, in QED, a = 0.2 or 0.3. 
QCD perturbation theory will fail at distances smaller than 0.2 fm, but not often 
and not by much; one should worry about one's simulation if it disagrees with 
perturbation theory at such distances. 

3.2 Keeping Up 

Improved lattice actions and operators can be quite complicated, and lately 
they have been changing rapidly. Unfortunately perturbative calculations of 
coupling constants and renormalization constants must be redone if any of the 
actions or operators involved is changed. Also the vertices and propagators 
can be quite complicated in a highly improved theory, making even one-loop 
calculations a chore. The perturbative work on NRQCD, which has a relatively 
complicated action, demonstrates that such calculations are possible [J10| . We 
need to systematize such calculations so that one-loop calculations, at least, 
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Figure 3: Clover coefficient c sw for SW quarks from Monte Carlo simulations, 
and from first-order perturbation theory using the plaqucttc mean link u p and 
the Landau-gauge mean link uq. Relative errors shown for the perturbation the- 
ory are ±la v . The lines indicate the possible range of the nonperturbative c sw 
(due to finite-a errors) around the parameterization given in (j| . 

can be completed in a short time for any new quark or gluon action. Otherwise 
we will find ourselves tied to older, probably much less efficient actions and 
operators by the lack of perturbative calculations for the new ones. Substantial 
progress in this direction is being made now, but more effort is needed. 

One approach that deals with this problem is to use Monte Carlo simula- 
tions to do the lattice perturbation theory [|TT]. On-shell quark form factors, 
renormalization constants, masses, etc. can be computed using simulations at 
very small lattice spacings, for example a « 0.002 fm at f3 = 9, and matched 
directly with continuum calculations. Since confinement is irrelevant at such 
short distances, coupling constants determined in such a matching can be fit by 
polynomials in ay (which is also measured in the simulation). In this way we 
obtain perturbative expansions for the couplings that can be used in simulations 
at more normal lattice spacings. 

4 Recent Developments 

4.1 Landau-Gauge Mean Link 

It is standard practice to approximate the Landau-gauge mean link uq by the 
plaquette mean link u p . When the lattice spacing is small, either of these is 
a reasonable choice. The two definitions differ, however, by almost 10% when 
a = 0.4fm (probably due mostly to perturbative differences in 0{a v )). While 
coupling constants obtained from the two definitions agree when perturbative 
corrections are included to all orders, tree-level or one-loop results can be quite 
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different when a is large. Three pieces of evidence suggest that it is better 
to stay with the Landau-gauge mean link uo- First the tree-level, tadpole- 
improved gluon action in Eq. (Q) shows 2-3 times less violation of rotation 
invariance (in the static potential) when uq is used rather than u p . Second, 
scaling violations in the charmonium groundstate hyperfine splitting are 3 times 
smaller when the tree-level a ■ B interaction is tadpole-improved with uq rather 
than Up Finally, the nonperturbative determination of the clover coeffi- 
cient c sw for SW quarks agrees more closely with tadpole-improved one-loop 
perturbation theory when uo is used, as is evident in Fig. ||. The same has also 
been observed at larger lattice spacings, with improved gluon actions, and also 
in improved quark actions like the D234c action discussed below Q . These 
examples suggest that uq might be the better choice generally, at least for large 
lattice spacings, if only tree-level or first-order corrections are included; there is 
little difference if corrections through second order are included. 



4.2 Anisotropic Lattices 

It is widely felt that coarse lattices are noisier than fine lattices. This is cor- 
rect, especially for correlation functions where signal/noise vanishes expontially 
quickly with increasing separation t. It is true despite the fact that the ratio of 
signal to noise is approximately independent of lattice spacing when t is large; 
noise overcomes signal at some fixed physical distance, independent of the lattice 
spacing. The problem on coarse lattices is that this fixed distance might be only 
one or two lattice spacings long, since the lattice spacing is large. Consequently, 
without large Monte Carlo data sets, the statistically usable part of a mass 
plateau, for example, may be too short to allow identification or verification of 
the plateau. A lattice with small lattice spacing allows us to identify the plateau 
at smaller distances, with less noise, but at increased cost from the small lat- 
tice spacing. Since signal/noise decreases exponentially with increasing t, while 
simulation cost grows only as a (large) power of lattice spacing, one might de- 
cide that smaller lattice spacings are ultimately more efficient. Much of the 
advantage of the larger lattice spacing is regained by careful design of sources 
and sinks, but a far more efficient approach is to reduce the lattice spacing at 
in the time direction while retaining a large spacing a s in spatial directions. 
Then more values of a correlation function can be measured at small t's, where 
relatively few Monte Carlo measurements are needed to suppress the noise. The 
small increase in simulation cost that results from a smaller a t is more than 
offset by the exponential improvement in signal-to-noise at the smaller t's. 

Anisotropic actions are easy to design. The standard improved gluon action 
becomes 15,0 



s n = -,3 £ (is) 

x,s>s' 
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12 u 6 s 12 
^ a s _-, f 4 P s . t 1 i? s 



(16) 



where m s and Ut are the (Landau-gauge) spatial and temporal mean links, re- 
spectively, and rj accounts for renormalization of the anisotropy by quantum 
effects. Since a t <^a s , we have not bothered to correct for order-a 2 errors in 
this action; thus the action extends at most one step in the time direction and 
has no gluon ghosts. An a 2 -improved version of Landau gauge is obtained by 
maximizing 



E 



-^ReTr\u il (x)--^—U ti (x)U tl (x + fi)\. (17) 

The anisotropic gluon action seems to have six parameters, but these appear 
in only three independent combinations, 

/? ( rja t \ u t 

u s , (18) 



u^ut ' V a s I u s 



which control the spatial lattice spacing, the anisotropy, and the size of the 
order a 2 corrections, respectively. These can be tuned numerically by choosing 
[3 and r)at/a s , then determining u t and u s self-consistently, and finally measur- 
ing rj. Alternatively one might compute the last three parameters using pertur- 
bation theory. In |l5|] simulations with very small lattice spacing are combined 
with first-order perturbation theory to obtain (prelimary) perturbative expres- 
sions for the case r\atja s = 0.5: 

u s = l-0.517//3+.17(5)//3 2 (19) 
u t = 1 - 0.105//? - .02(l)//3 2 ( 20) 

?7 = 1-0.114(7)//?. (21) 

These expressions give good results even at a = .4fm, as shown in Table 0. 

In Fig. ^ we show the static potential computed on an anisotropic lattice from 
Wilson loops involving temporal and spatial links, and from loops involving just 
spatial links. Shifted to account for renormalization, the two potentials agree 
well for all r, as they should. 

Anisotropic lattices with a t a s are useful for any simulation involving high- 
energy states, including glueballs, excited hadronic states, and high-momentum 
final states in form factors. The dominant lattice errors in relativistic heavy- 
quark physics are powers of a t Mq, where Mq is the mass of the heavy quark; 
these are greatly reduced on highly anisotropic lattices, which, therefore, pro- 
vide an alternative to the NRQCD and Fermilab approaches O] . A variety of 
quark and gluon actions designed for anisotropic lattices are currently under 
development. 
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(3=1.8 


(3=1.34 


U s 


P.Th. 


.77(2) 


.71(3) 




M.C. 


.75 


.69 




P.Th. 


.94(1) 


.91(1) 




M.C. 


.93 


.91 


V 


P.Th. 


.94(1) 


.91(2) 




M.C. 


.94(1) 


.90(1) 



Table 1: Monte Carlo and perturbation theory values for parameters in the 
anisotropic gluon action with rjat/a s = 0.5. Results are given for two lattice 
spacings (approximately 0.25 fm and .4fm), and are preliminary. 
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Figure 4: Static potential computed on an r ^!nfeotropic lattice with rja t /a s = 0.5 
and (3 = 1.8. Potentials obtained from spatial-temporal loops and from purely 
spatial loops are shown. Results are preliminary. 
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4.3 Improved Gluons 

The plaquette-plus-rectangle gluon actions, Eqs. (|l|) and (|l7|), have been used 
by several groups in recent years for simulations of the static potential, QCD 
thermodynamics, charmonium, the light-hadron spectrum, and so on. In every 
application area, except perhaps the light-hadron spectrum, the improved action 
has led to significantly improved results. 

One of the most striking recent applications is in jl6| , a high-statistics study 
of several different glueball masses at a wide range of lattice spacings. What is 
striking is that this thorough study was carried out entirely on desktop worksta- 
tions, an achievement possible only through the use of coarse, highly anisotropic 
(5:1) lattices. All the glueball masses extrapolate to the values determined ear- 
lier with more traditional methods; and all but one show relatively little de- 
pendence on the lattice spacing, less than 10%, out to lattice spacings of order 
0.4 fm. Only the ++ glueball shows strong a dependence; its mass has a 20- 
25% dip around 0.3 fm. These are the largest measured finite-a errors in gluonic 
lattice QCD, and therefore the ++ glueball is of particular interest to designers 
of improved gluon actions. 

A dip in the ++ mass is also seen in SU(2) QCD |l7|], although it is smaller, 
only about 14%. The authors of the SU(2) analysis noticed that the depth of 
the dip is reduced almost by half when the Landau-gauge mean link uq is used, 
in place of the plaquette mean link u p , in the tadpole-improved tree- level gluon 
action. This result is consistent with the observation that the static-quark 
potential is more rotationally invariant with uo; both argue that the correct 
order a 2 correction in the action is stronger than that obtained using u p . The 
SU(3) study used u p ; it should be redone with uq instead. Better still, we should 
calculate the full 0(a 2 av) corrections to the anisotropic gluon action; work on 
this has begun. A small glueball radius would explain the strong a dependence 
of the ++ mass; another explanation is that the dip in the mass is caused by a 
phase transition nearby in the fundamental-adjoint coupling constant plane of 
the gluon action. 

4.4 Improved Wilson Quarks 

Many of the talks at this meeting showed results for the SW quark action, 
which is an order a improved version of the original Wilson lattice action. Re- 
sults are clearly superior when the SW action is used instead of the Wilson 
action, provided the clover coefficient is cither tadpole improved or tuned non- 
perturbatively; and generally SW and Wilson simulations extrapolate to the 
same answer. Based on this work, the SW action should be accurate to within a 
few percent for light-quark physics out to lattice spacings of order 0.15-0.2 fm. 

Larger lattice spacings require a more accurate action. A simple construction 
that yields such an action begins with the continuum lagrangian: 

C = i> c M c ^ Cl M c = D • 7 + m c . (22) 
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Still working in the continuum, we introduce a field transformation, 



ip c = n^p (23) 

fl 2 = 1 - \ra(D ■ 7 - m c ) (24) 
to obtain a new lagrangian ipM^ where 

m = riM c n 

= D ■ 7 + m c (l + \ram c ) (25) 
-\raip 2 + \ga ■ F) . (26) 

The third term in M , introduced by the field transformation, breaks the dou- 
bling symmetry when we discretize the action. Allowing errors only of order a 4 
and higher, we obtain the D234c discretization of M [H : 

a 2 

M D234 ee A • 7 - c 3 — A< 3 > • 7 + m (27) 

2 V 2 12 

where m = m c (l + iram c ), A!,™" 1 is the lattice discretization of the n'th gauge 
covariant derivative D™, and F^ v is an a 2 -improved clover discretization of the 
gluon field strength. All operators are tadpole improved; at tree level the cou- 
plings Cp = Cs = C4 = 1. 

We have tested this action on lattices with a = 0.25 fm and 0.4 fm JlJ]. We 
limited our study to quark masses near the strange quark mass. Thus we avoid 
the inessential complications of chiral extrapolation. Hadrons made of quarks 
with larger masses are generally smaller and heavier; consequently discretization 
errors are probably larger for the strange quark than for the other light quarks. 

The leading correction in Md234c beyond those in the SW formalism is 

-c 3 yE A ! 3 V (28) 

This violates Lorentz invariance; it cancels similar violations in the leading-order 
terms. Thus we tested this correction by computing 

C '(p) S £ '- 1P) p ' £ ' (0) (») 

where Eh(jp) is the energy of hadron h at three-momentum p. In a Lorentz in- 
variant theory, c(p) should equal 1, the speed of light, for all p. Results at 0.4 fm 
are shown in Fig. || for both the D234c and SW actions. D234c is dramatically 
superior, deviating from 1 by only 3-4% at p = 0. Also the dispersion relation 
is accurate to within 10% out to three- momenta of order 1.5 /a. 
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Figure 5: Speed of light squared from the quark dispersion relation using the 
D234c and SW actions, for pseudoscalar ("7r s ") and vector ("0") mesons made 
of strange quarks. Simulations are on isotropic lattices with a«.4fm and an 
improved gluon action. Fits shown are to c + C4(pa) 4 for D234c and c + c 2 (pa) 2 
for SW. 

We could tune the coupling C3 in the D234c action nonperturbatively until 
c 2 (0) = 1, thereby restoring Lorentz symmetry through order a 2 . In the tadpole- 
improved theory at 0.4 fm, however, c 3 = 1 is so close to the correct value there is 
no point. There is no value of C3 for which both the pseudoscalar and vector have 
exactly c 2 = 1 ; uncorrected a 4 errors are larger than the radiative corrections to 
C3 at this lattice spacing. 

There are two different definitions of a hadron's mass in lattice simulations. 
One is the "static mass," Eh(0), and the other the "kinetic mass," which equals 
E h {0)/c 2 (0). In D234c these two definitions agree to within 3-4%. In SW they 
differ by 40-60%. The kinetic and static masses in either formalism must be 
equal for zero mass mesons, because of the axis-interchange symmetry of the 
actions. The deviations seen here are because the strange quark is relatively 
massive at 0.4 fm; the <f> mass, for example, is 2.1/a with D234c. 

Another test of the D234c action is to vary the parameter r. The r-dependent 
terms in the D234c action are approximately redundant at low momenta, since 
they originate from a field transformation. We ran simulations both at the 
standard value r = 1 and at r — 2/3 to check this. With D234c the hadron 
masses changed by no more than 1%. Reducing r in this way lowers one of the 
ghost branches of the quark dispersion relation from an energy near 2/ a to one 
near 1/a. The lowering of the ghost changed little in the spectrum; our results 
are not strongly affected by ghost artefacts. 

With SW, where the r terms are less highly corrected and therefore less 
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redundant, we found shifts of order 1-5% when we varied r as above. 

The coefficients cf, C3, and C4 are all renormalized by residual quantum 
effects, even after tadpole improvement. We varied these couplings to assess 
the importance of the potential corrections. Only cf had a strong effect on the 
static hadron masses. D234c with cp = 1 gives scaling errors in the <j> mass of 
order 5% at 0.25 fm and 15% at 0.4 fm; but a radiative correction of 0.5ay in 
Cp, like that in the SW action (Eq. (|l4|)), would shift these masses by 4-8%. 
Thus to achieve precisions of order a few percent at such large lattice spacings, 
we require a calculation of the order ay correction to c sw . This calculation will 
be completed in the very near future. 

The highly corrected D234c action is far superior to SW for states with 
high energies, either because p is large or because a is large. This makes it 
potentially useful for simulations with very large lattice spacings, large quark 
masses, or high-momentum states (for example, high-energy weak interaction 
form factors). Anisotropic versions of D234c have also been tested (lj, |l). 

4.5 Improved Staggered Quarks 

The standard staggered-quark action, like the SW action, has leading errors of 
order a 2 . Some attempt has been made to correct the a 2 errors jl8| , but this 
effort included only the obvious correction to the derivative in the action. Sur- 
prisingly, there are, in addition, four-quark interactions that must be included 
at tree-level in order a 2 . Tree-level improvement of a lattice action requires 
that we match on-shell, tree-level, lattice scattering amplitudes with continuum 
QCD at low external momentum. Normally this implies that the lattice theory 
has no tree-level four-quark interactions since continuum QCD has none: in a 
tree-level quark-quark scattering amplitude, low-momentum external quarks im- 
ply a low-momentum internal gluon, and therefore these amplitudes will match 
automatically if the quark-gluon vertices match. This reasoning is incorrect 
for staggered quark formalisms. For example, a pair of low-momentum quarks 
can scatter exchanging a high-energy gluon with momentum q — (n/a, 0, 0, 0), 
creating a pair of low-momentum quarks of a different flavor. Such flavor chang- 
ing quark-quark amplitudes are absent in QCD and so must be eliminated by 
adding correction terms to the action. Since the gluon is highly virtual in a 
flavor-changing process, the correction term is a four-quark operator with a co- 
efficient of a constant times a 2 av(q g ), where the constant is computed from the 
scattering amplitude and q g is the gluon momentum. 

The four-quark corrections obviously will greatly reduce the flavor-symmetry 
violations observed in hadronic spectra from staggered quarks, since these are 
proportional to a 2 . Supporting evidence also comes from the numerical ex- 
periments described in |Q. In these experiments, the link operator in the 
standard staggered-quark action was replaced by a smeared link. Smearing the 
links inhibits the exchange of high-momentum gluons by lowering the effective 
ultraviolet cutoff of the gluons; but the flavor-changing interaction requires ex- 
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changed gluons with the highest momenta on the lattice, ±7r/a in one or more 
directions. Thus smearing should dramatically reduce the flavor splittings; and 
it does . 

If the four-quark operators are computed and can be included in a practical 
action, the staggered-quark formalism will be more accurate than the SW ac- 
tion. Adding the four-quark corrections may also significantly reduce the size of 
ay corrections to the mass renormalization for staggered quarks. This correction 
term has normal size for r~l, but becomes quite large in the staggered-quark 
limit, r = 0. Thus the large size in this limit seems to be related to flavor mixing. 
A reduction of this renormalization is important for the determination of the 
strange-quark mass. 

5 Conclusions 

Techniques for simulating on coarser lattices are rapidly gaining acceptance 
within the lattice QCD community, and work at lattice spacings in the range of 
0.1-0.2 fm is becoming widespread, again. It is likely that some combination of 
perturbatively improved actions, nonperturbative improvement, perfect actions, 
anisotropic lattices, Monte Carlo Renormalization Group, and the like will allow 
us to simulate accurately and efficiently on still coarser lattices in the very near 
future. 
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